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We study a spin 1 Bose condensate small enough to be treated as a single magnetic 'domain': a 
system that we term a microcondensate. Because all particles occupy a single spatial mode, this 
quantum many body system has a well defined classical limit consisting of three degrees of freedom, 
corresponding to the three macroscopically occupied spin states. We study both the classical limit 
and its quantization, finding an integrable system in both cases. Depending on the sign of the 
ratio of the spin interaction energy and the quadratic Zeeman energy, the classical limit displays 
either a separartrix in phase space, or Hamiltonian monodromy corresponding to non- trivial phase 
space topology. We discuss the quantum signatures of these classical phenomena using semiclassical 
quantization as well as an exact solution using the Bet he ansatz. 



I. INTRODUCTION 

A spin 1 Bose condensate can be regarded as an un- 
usual kind of magnet. As with a magnet, we can con- 
sider a small system that constitutes a single magnetic 
'domain'. The study of the dynamics and excitation spec- 
trum of a spin 1 condensate in this regime - which we 
term a microcondensate - is the subject of this paper. 

A number of recent experiments have demonstrated 
the relevance of this 'single mode' approximation in ex- 
periments using condensates of ^^Rb or ^^Na atoms [T]-[3]. 
This regime occurs when the size of the condensate is less 
than the spin healing length over which the spin state of 
the particles can change significantly. This allows us to 
consider the spin dynamics in a single spatial mode. In 
terms of the operators Aj^, m = —1,0,1 that cre- 
ate and destroy particles in the three spin states of that 
mode, the Hamiltonian for a system of N particles is 

HsMA = ^:N':+^:S-S:+H^. (1) 

Here the colons denote normal ordering and TV and S are 
respectively the operators of total number and spin 

1 

N= ^ AI^A^ S= ^ Al^Smm'^rn'^ (^) 

m= — l m^m' 

where S^n are the spin 1 matrices. The energies eo and 
62 parametrize the strengths of the density-density and 
spin-spin interactions, the most general situation for spin 
1 [4^ . i^z describes the Zeeman energy, and includes both 
linear and quadratic contributions 

1 

Hz= ^ [Pm + qm"] A^ (3) 

m— — l 

The investigation of the Hamiltonian Eq. ([T]) is the main 
goal of this paper. It is integrable both at the classical 
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and quantum level. We shall show that this rather simple 
system displays a variety of interesting properties that 
merit further experimental study. 

When the number N is large we expect on general 
grounds that the quantities A^, A\^ can be treated as 
classical amplitudes A^, A^, with a classical Hamilto- 
nian given by Eq. ([T]), and Poisson brackets that reflect 
the more familiar quantum commutators 

{^m^^m' }pB ~ J^^mm'- (4) 

The resulting equations of motion are nothing but the 
'zero-dimensional' Gross-Pitaevskii equation. This me- 
chanical system has three degrees of freedom (so that 
the phase space is six dimensional, corresponding to the 
three complex quantities A^, m = —1,0, +1), and three 
commuting conserved quantities: N ^ Sz-, and the energy 
HsMA itself (note that is not conserved because of 
the quadratic Zeeman effect). As a result, the system is 
integrable in the sense of Liouville, a notion that will be 
reviewed in Section BlI Bl 

We will see that the character of the classical motion 
depends essentially on the sign of q = qje^- Since q is 
typically positive (though an effective negative q can be 
induced by microwave dressing [5l |6]), this means that 
g > is most naturally realized for antiferromagnetic in- 
teractions e.g. in ^^Na and g < in the ferromagnetic 
case e.g. ^^Rb. For a graphical illustration of the differ- 
ence between these two cases, the reader should compare 
the spectra shown in Figs, [s] and [s] (note that the spectra 
shown are all plotted with 62 = 1, so the spectrum should 
be inverted in the ferromagnetic case). The way in which 
the classical motion influences the quantum spectrum is 
one of the main themes of this paper. 

The structure of the remainder of this paper is as fol- 
lows. In Section |ll| we recast the Hamiltonian Eq. ([T]) 
in terms of hyperbolic (or 6'[/(l,l)) spins, which pro- 
vides a convenient framework for the analysis of both 
the classical and quantum systems. Section [lll| presents 
a detailed analysis of the classical problem and its semi- 
classical quantization. After describing the features of 
the classical dynamics in qualitative terms and explain- 
ing the relation to the mean-field phase diagram, we in- 
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troduce some ideas from the theory of classical integrable 
systems, notably action- angle coordinates, that are then 
applied to the system of interest. It is in this section that 
we meet the phenomenon of Hamiltonian monodromy^ a 
topological obstruction to t he glo bal existence of action- 
angle coordinates (Section |lIID[ ) [7H9]. This material 
may be unfamiliar to many readers, so we have tried to 
be pedagogical in our presentation (other introductions 
suitable for physicists may be found in the appendices 
to Ref. [10], see also Ref. [JlJ). Finally Section IV is de- 
voted to a very different approach, in which the quantum 
Hamiltonian is solved directly using the Bethe ansatz. 
The Bethe ansatz equations were given in Ref. [12 , based 
on mapping to a type of Gaudin model written in terms 
of the hyperbolic spins. Here we re-derive these equations 
by a different method and describe their solution, with a 
particular emphasis on the properties already uncovered 
in Section [TTTl 

We close this introductory section with a couple of 
comments on some related work. The low energy physics 
of spin 1 condensates with antiferromagnetic interactions 



(also called 'polar' condensates) has often been studied 
in terms of an effective 'rotor' description, akin to the low 
energy description of the Neel state of an antiferromagnet 
[13^. Refs. [iT, 15^ recently extended this description to 
the full spectrum in the single mode approximation. The 
rotor formulation is quite different from the approach 
pursued in this work, however. Finally, a model identical 
to Eq. ([T]) was studied numerically, and the existence of 
monodromy pointed out, in Ref. [16] in a very different 
context. Our goal here is to provide an analytic descrip- 
tion. 



II. FORMULATION IN TERMS OF 
HYPERBOLIC SPINS 

Let us write Eq. ([T]) explicitly in terms of the boson 
field operators. For clarity we drop the terms propor- 
tional to eo (density-density interaction) and p (linear 
Zeeman effect) as these couple to conserved quantities, 
giving the reduced Hamiltonian 



red 
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2N 



Si + 2 (a\aIMo? + {AlfA,A_,) + 2 (^4^0 - (a\A, + aI,A 



'a\a_^+a1^a_ 



(5) 



In this form the Hamiltonian is still somewhat indi- 
gestible. A considerable simplification was achieved in 
Ref. [T2j by introducing the variables 



K+ 
Bo 
B. 



1 

2 

a'' a'' 

^1^-1' 



a\a^ + a1^a_-^ + 1 

K_ = A,A 



-1' 



0^0 



1 

4' 



(6) 



-i(44),i?_ = -i(^Ao 



We make two important comments about this choice of 
variables. Firstly, they are invariant under the rotations 
about the z-axis generated by the conserved quantity Sz 
(note that from now on we will we drop the tildes on 
N and Sz, trusting that no confusion between the quan- 
tum operators and their eigenvalues will result), and are 
therefore suited to exploiting this symmetry of the sys- 
tem. Secondly, they constitute two representations of the 
non-compact group SU{1^ 1), obeying the relations 



[Ko,K±]=±K± 
[K+,K-] = -2Ko 



(7a) 
(7b) 



and likewise for the {Bq^ B-^, B-} variables. The differ- 
ence from the more familiar SU{2) algebra lies in the 
minus sign in Eq. (|7b[), leading to the quadratic Casimir 



operators 

Ck = Kl-\ {K.K+ + K+K.) = \ [Si - 1) 



Cb=BI-\{B_B++B+B.) 



3 

16' 



(8) 



In the first case the representation is fixed by specifying 
the value of Sz = A\A-^ — A^_-^A_-^. In terms of the 
occupation number basis \Ni)-^ of the m = ±1 

states the highest weight state is 

where the Bargmann index = ^ (l*^'^! + 1) is the eigen- 
value of Kq for the highest weight state. Note that 
Ck = i^K {i^K — !)• Repeated application of gen- 
erates the states 



|n, uk) k = \/(2z^K + n)(n + 1) \n + 1, vk) k 
, . ^\\S.^n),\n)_,, iiSz>^ (10) 

with Kq eigenvalue 

Kq |n, vk)k = (^ + ^k) \n, vk) • 
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For the {Bq, B-} variables ('one mode representa- 
tion') things are shghtly different. There are just two 
representations, with highest weight state |0)q or i.e. 



either no particle or one particle in the m 
These have Bq = so 



\0,iyB 

10,1^3 



1/4)b 
3/4)5 



|0>o 
|l)o 



both giving Cb = {^b — 1) = 
determined by the parity of N — Sz 



16' 



state. 



(11) 



This index is 



J^B 




for AT 



Sz even 
Sz odd 



(12) 



Repeated application of B-^. generates the states 
B+ \n, 1^3)3 = Vi'^'^B + n){n + 1) \n + 1,vb)b 



\n,VB)i 



_ j(-l)"|2n)o, 



if vb 



(13) 



(-l)"|2n + l)o, ifi^B = 



In the classical limit we can think of the Casimirs of 
Eq. ([s]) as specifying a hyperboloid in the space of values 
{K,,Ky,K,) with [17. 



K-) 
-K-) 



(14) 



and similarly for the B representation. Thus we will refer 
to these variables as 'hyperbolic spins'. In this language 
the Hamiltonian Eq. ([5| takes the form 
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[Sl + SBoKo - 4B_K+ - 4B+K_] 



2qKo. 
(15) 

We have dropped a term — 62(1+2^)— ^, constant by 
number conservation. In Eq. (IT5| the conservation of 
particle number N = 2{Ko + Bo)~- | is also manifest as 
a symmetry under rotations about the 2:-axis in the space 
of the hyperbolic spins (see Eq. ([l4|). We readily see that 
in terms of the original Aj^ degrees of freedom, this is just 
the overall phase conjugate to the total particle number. 



We emphasize that Eq. ([15| is identical to the origi- 
nal Hamiltonian Eq. ([5|. The symmetry under rotations 
about the z-axis (equivalently, conservation of Sz) allows 
us to write the original problem with three degrees of 
freedom as a problem with only two degrees of freedom. 
The price we pay is that the relevant variables are the 
less familiar hyperbolic spins. 

This formulation of the problem turns out to be a very 
convenient starting point for the analysis of both the clas- 
sical and quantum dynamics of the spin 1 condensate. 
III. THE (SEMI-)CLASSICAL LIMIT: 
REDUCTION AND MONODROMY 

A. Qualitative features of reduced dynamics 

In this section we analyze the classical mechanics of 
the system described by Eq. ( 15 ). On general grounds we 



expect this to be a good description when the number of 
particles N is large, in which case the operators A^, Al^ 
can be treated as classical amplitudes A^, of 'order 



Thus the B and K variables are 0{N), as is 5'^, 
so in this limit we take 



Ck 



4 
~ 0. 

iV ~ 2 (Bo 



(16) 



and the 0(1) terms in Eq. ([6| can be dropped. Thus the 
B degrees of freedom are restricted to a cone Bq = 
Note that Ko,Bo< f. 

We can exploit the symmetry under rotations about 
the z-axis in hyperbolic spin space, corresponding to con- 
servation of particle number, by taking 5+ to be real and 
positive. Then we have 5+ = Bq = N/2 — Kq and we 



can eliminate the B degrees of freedom from Eq. (15) to 
give the classical Hamiltonian 



-f^class — 
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2N 



[Sl^8{N/2-Ko){Ko-K, 



2qKo. 
(17) 

If the reader finds these manipulations too cavalier, an 
alternative is to use conservation of N to write the quan- 
tum Hamiltonian Eq. (15) as 
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27V 



St ^8 



N 
~2 



2^{N - 2Ko){N - 2Ko + 1)K+ - 2K_^{N - 2Kq){N - 2Ko + 1)1 + 2qKo, (18) 



which reduces to Eq. (17) when the 0(1) terms are Now of course, the semi-classical limit is the domain 
dropped. of validity of the conceptually simpler framework of the 
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Gross-Pitaevskii equation (see Ref. [18], for instance). 
Why go to all this effort? The answer is that we have 
exploited the conservation laws of the problem to reduce 
the dynamics to a single degree of freedom, which allows 
us to say a good deal about the character of the motion 
without detailed calculations. We will refer to the hyper- 
boloid of fixed Sz with Kq < N/2 as the reduced phase 
space. 

For the remainder of this section we work with 'per 
particle' quantities in lowercase: thus kx = K^/N^ 
(ai,ao,a_i) is the condensate wavefunction normalized 
to unity, the dimensionless energy per particle is = 



^ciass/^e2, and the dimensionless quadratic Zeeman 
shift is g = g/e2. Thus 



h=]^sl^ 2(1 - 2/co) (^0 - kx) + 2qko (19a) 



= ±2x1 kl 



(19b) 



Let us begin by discussing the level sets of the Hamilto- 
nian function Eq. (19a). This is just a quadratic form, 
which we may write 



h=\sl^2 (fco - 1/2, K + D) (f ;) {,'1 1)) + i 



(20) 




► 



FIG. 1. Example of the geometric construction for = 0.5 
giving the trajectory (dark line) as the intersection of a hy- 
perboloid, the reduced phase space of fixed Sz = 0.2, and the 
surface of fixed energy h — 0.7. 



The surface defined by = const, is a hyperbolic sheet 
whose normal lies in the {ko^kx) plane. The asymp- 
totes of this surface are the planes ko = ^ and ko = 
kx -\- ^. Theses planes correspond to the vanishing of 
the quadratic form in Eq. (20) and hence to the energy 
h 



2^? - 



-q. The trajectories of the reduced Hamiltonian 
are given by the intersection of the hyperboloid of fixed 
Sz with the surface of fixed h (see Fig. [T]). 

The character of the trajectories depends critically on 
the sign of and we discuss the two possibilities in turn. 



^ > 



kx -\- ^ slices 
Thus, the in- 



For < q < 2 the asymptotic plane ko - 
through the hyperboloid below ko = 1/2. 
tersection of the asymptotes with the hyperboloid forms 
a separatrix at energy h = ^s^ -\- q (see Fig. [5]). Since 




0.0 



FIG. 2. For positive q (here q = 0.5) the asymptotic planes 
/co = |, A^o = A:^ + f (black dotted line) of the hyperbolic 
sheets of constant energy cut the hyperboloid of fixed Sz — 
0.2 (blue dashed line), forming a separatrix between energies 
greater (red line) or less (green line) than h — \s1 -\- q. 



6o = ^ — /co it is clear that the separatrix corresponds to 
a trajectory for the {hx^ by^bz) variables that passes over 
the tip of the cone bo = |&+|. 

For a fixed 5^, only a certain range of values of h are 
allowed, corresponding to non-zero intersection of the 
h = const, surface with the hyperboloid. The endpoints 
of this range correspond to the stationary points of h on 
the hyperboloid of fixed Sz^ which are obtained from the 
equations 



d{h-Xsl) 

dko 
d{h-Xsl) 

dkr 



4:k^ 



2 - 8/co + 2(7 - SXko = 



(21) 



2(1 - 2ko)±SXkx = 0, 



where A is a Lagrange multiplier. This pair of linear 
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FIG. 3. Spectrum of Eq. ([5| for N = 30 particles with q = 
0.3. The sohd hnes correspond to the boundaries given by 
{sz{kx,ko),h{kx,ko)) evaluated on kx,o{X) from Eq. (22) for 
the ranges in Eq. ( |23| . The dashed line corresponds to the 
separatrix /i = ^sl -\- q. 



FIG. 4. (For negative q (here q = —0.5) there are two distinct 
classes of trajectories, that either do (green line h < 0) or do 
not (red line h > 0) encircle the apex of the cone correspond- 
ing to Sz = 0. 



equations gives ko as a function of A 



l + 2A(l + (7) 

2(1 + 2A)2 V l + 2A-(7. 



(22) 



These values can be inserted into h{kx^ko) and 
Sz{kx^ ko) = 2y^/cq — /c^ to yield a parametric curve. 

There are two ranges of A that correspond respectively 
to the upper and lower limiting values of h for each Sz. 

q-2 ^ q-1 

upper range: < A < 

^ 4 + 2^ 2 (23) 

lower range: < A < oo 

For q < 1 this does not include the full lower boundary 
because for < A < the solution in Eq. ([22|) corre- 



sponds to values ko > 1/2. In this case part of the lower 
boundary is given by the separatrix (see Fig. [3| . Once 
q > 2 the upper boundary is given by the separatrix. 



2. q<0 



In this case the asymptotic plane ko = + f ^^^^^ cuts 
the hyperboloid. The lower energy boundary corresponds 
to trajectories on the asymptotic plane ko = ^ with h = 
-\- q. The upper boundary is found as in the ^ > 
case. 

An interesting situation does arise, however, for Sz = 0. 
Here the hyperboloid of fixed Sz becomes a cone ko = 1 
and we can distinguish trajectories based upon whether 
or not they encircle its tip (see Fig. [4]). From Eq. (19a) 
we see that the critical value h = separates these two 
regimes, with h < encircling the tip and h > not. 
The feature at the origin of the spectrum shown in Fig. |5] 




-30 -20 



FIG. 5. Spectrum of Eq. ^ for N = 30 particles with 
q = —0.5. The solid (upper boundary) is found as in the ^ > 
case, and the lower boundary (dashed) corresponds to trajec- 
tories on the asymptotic plane /co = | of energy /i = -\- q. 
The feature at the origin is a signature of the phenomenon of 
monodromy. 



is a signature of the phenomenon of monodromy to be 
discussed shortly, which is related to this distinction. 

When q < —2 the origin passes through the up- 
per boundary (note the behavior of the lower limit in 
Eq. J23|), after which the upper boundary is cusped (see 
Fig. |6|. The relation of these changes in the morphology 
of the spectrum to the mean-field phase diagram will be 
discussed in the next section. 
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FIG. 6. Spectrum of Eq. ([5| for AT = 30 particles with 
q = —2.2. For q < —2 the origin passes through the upper 
boundary, which becomes cusped. 




q 



FIG. 8. The ground state phase diagram for 62 > 0. The 
shaded region is where the separatrix (dashed hne in in Fig.|3| 
is the lowest energy. 
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FIG. 7. The ground state phase diagram for 62 < 0. The 
vectors denote the spinor (ai,ao,a_i). The shaded region 
corresponds to the dark upper boundary in Fig. |5] and Fig.|6] 



S. Relation to mean-field phase diagram 

The mean-field phase diagram for the ground state in 
terms of the linear and quadratic Zeeman shifts p and q 
was given in Ref. [19 . In fact, p is more properly regarded 
as a Lagrange multiplier used to find the ground state for 
fixed Some of the features discussed in the preceding 
two sections can be related to the structure of the phase 
diagram. 

Let us begin with the ferromagnetic case, correspond- 
ing to 62 < 0. The phase diagram is shown in Fig. [7[ 
Positive q (the most physically relevant case) corresponds 
to g < 0, and the ground state to the largest value of h 
(since i^ciass = h/e2N). Thus as p increases we move 
along the dark upper boundary in Fig. [5] and Fig. [6j cor- 



responding to the shaded region in Fig. [7[ The difference 
between q > —2 and ^ < — 2 is that in the latter case 
a finite p is required before we depart from = be- 
cause of the cusp in Fig. [6) The transition corresponding 
to the lower boundary of the shaded region in Fig. [7| is 
the zero-dimensional analog of the transition discussed 
in Ref. [20 . 

In the antiferromagnetic (or polar) case 62 > 0, we 
follow the lower boundary in Fig. [3] Because the cusp 
is always present for ^ > a finite p is always required 
to depart from Sz = 0. The negative curvature of the 
black lower boundary means that the ground state jumps 
straight to a point on the dashed (separatrix) line when 
\p\ exceeds the critical value ^/2qe2. For q> \^ the sys- 
tem jumps straight to full polarization at p = g + ^ (see 
Fig.[8|. 



B. Classical integrability and its consequences 

We began with a system with three degrees of free- 
dom and have reduced the problem to the motion along 
constant energy contours on a hyperboloid. This reduc- 
tion used the conservation of N and Sz^ which together 
with the energy constitute three independent commut- 
ing conserved quantities. On general grounds, once we 
fix the values of these conserved quantities, the system 
must move on a three dimensional submanifold of the 
six dimensional phase space. The Liouville-Arnol'd the- 
orem provides more detail, telling us that this submani- 
fold is in fact a three-torus = 5'^ x 5^ x S'^ [21]. The 
three circles correspond to the overall phase of the spinor 
(^A\^ A^^ A_\) (conjugate to N\ the angle of rotation 
about the z-axis (conjugate to Sz\ and the closed tra- 
jectories of constant energy on the reduced phase space 

The constructive part of the Liouville-Arnol'd provides 
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a distinguished set of coordinates for this torus, known 
as action- angle coordinates. As this idea plays an im- 
portant role in what follows, we explain it in some de- 
tail. Each of the conserved quantites (we denote them 
by i = 1,...A^, for degrees of freedom) gener- 
ates a 'time' evolution on the torus, by using each in 
place of the Hamiltonian = Fi in Hamilton's equa- 
tions. Since these flows commute, the trajectory of a 
point X on the torus under their combined action can 
be written . . . ,tAr), with the first argument corre- 
sponding to the usual time evolution. Now the evolution 
in each variable is not in general periodic, but rather 
quasiperiodic, consisting of N incommensurate frequen- 
cies. The set of values of t = (ti,...,tAr) for which 
x{ti, . . . , t^) = ^(0, . . . , 0) is a lattice (the period lattice) 
consisting of integer linear combinations of some basis 
set ei,...,e7v, . We can then define the angular vari- 
ables in terms of the reciprocal lattice vectors €i, . . . , eat 
satisfying • ej = 27TSij 



ei • t, 



(24) 



which increase by 27r as we advance one unit in each of the 
lattice directions. These provide a natural parametriza- 
tion of the torus. Their time evolution due to the Hamil- 
tonian is particularly simple: 



(€i)i, i = 



l,...,Ar. 



(25) 



Note that for a given period lattice, the lattice vectors 
ei, . . . , e^v are not unique, resulting in an arbitrariness 
in the angles that will be important in the following. 
The final step is the introduction of the actions 



1 

2^ 



(26) 



where p and q are the canonical momentum and position 
variables, and the integral is taken around the i^^ circle of 
the torus. The definition Eq. (26) gives the li in terms of 



the conserved quantities Fi^ as specification of the latter 
fixes the torus. The relation may be inverted to give F^, 
notably the Hamiltonian, in terms of the 7^. 

The actions are conjugate variables to the angles in- 
troduced above, so that 



(27) 



Ik^Ii fixed 



The evolution generated by each of the conserved quan- 
tities corresponds to a matrix of angular 'velocities' 



/ IkT^h fixed 

or equivalently the period lattice vectors 

= (§-) 

V^^J / Fk^Fj fixed 



(28) 



(29) 



While one can add arbitrary constants to the actions 
without changing anything, there is more freedom in the 
choices of angles, where we may redefine 



(Pi + Ai(/i,...,/ 



N) 



(30) 



There is a close analogy to the gauge transformations 
familiar in quantum mechanics. In more mathematical 
terms, the phase space of an integrable system has the 
form of a fiber bundle, with the base manifold consisting 
of the space of conserved quantities F^, and the fibers be- 
ing the tori. As we move around the base manifold, we 
can change the definition of the angular variables on each 
torus arbitrarily. As in the Aharonov-Bohm effect, inter- 
esting things can happen when we move around a circuit 
containing a singular point, leading to multivaluedness 
of the angles. This provides one view on the monodromy 
that we will discuss in Section [IHDI 

From these generalities we return now to the system 
of interest. In this case the canonical form p • (iq be- 
comes I A'^dAjn — A^dA"^ (we set h = 1 for the 
remainder of this section). Now we choose the following 
parametrization for the spinor (for > 0) 



5:cosh-e-WWx) 
2 




(31) 



Note that x is the overah phase, while <p describes rota- 
tions about the z-axis generated by Sz, and 

K+=AIA*_^ = ye'*sinh6' 



\A_iY) = — cosh( 



4 ' 



so that 9 and xf) parametrize the hyperboloid. The logic 
behind this choice becomes apparent when we compute 
the canonical form 



^ Al^dAra - AradA*^ =Ndx + S, 
m 



cosh Odt/j, (32) 



where N is the total number of particles N = |^mP- 
In this way we can compute the three actions correspond- 
ing to the three circuits: x • ^ ^tt, (j) : ^ 27r, and the 
trajectories on the reduced phase space discussed in Sec- 
tion HI A Using Eq. ( 26 ) we find that the corresponding 



actions are 



h 
h 



N 

— ^ i 
2^ 2 / 



(33) 



cosh d^l) 
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The last formula has a geometrical interpretation in 
terms of the area of the hyperboloid enclosed by the 
trajectory. This is not the area induced by the usual 
Euclidean metric, but rather by the 'Minkowski' metric 



dKi - dKi - dKt = 



{de^ - sinh^ e di^"^) , (34) 



(c.f. the quadratic term in Eq. ( 15 )) with the correspond- 
ing area element 



dA 



sinh 9 dO dip 



Thus we have 



h = 



ttS, 



(35) 



(36) 



This interpretation will prove useful in the next section 
when we investigate the behavior of this action. 




C. Properties of the action 

Let us now investigate the properties of the action Is. 
Note that we are going to continue assuming > to 
avoid a rash of modulus signs. Du e to the symmetry of 
the Hamiltonian function h in Eq. (19a), it is convenient 



to work not with the coordinates and introduced 
in the preceding section, but instead with the half-plane 
model of hyperbolic geometry [22], illustrated in Fig. [o] 
The variables on the half plane are 



2 A^7i 



ko kx 
ko - kz 



(37) 



In terms of these variables fc, = 4f - while 

y z z 



ko 



knr — 



Sz 


z 






2" 

+ - 


4 






z _ 


Sz 


z 




y'\ 


2" 


4 




z^ J 


z 



(38) 



The equation relating y and z given h and is 

(5, - qzl2) = ^z^-z\2h/s,)^2z (2 + q)-4s, = Ps{z) 

so that 



y = ^\ 



Ps{z) 



{s, - qz/2) 



The area element is 
dA 



dydz 



(39) 



(40) 



FIG. 9. Trajectory on the reduced phase space in the half- 
plane model. The trajectory is first projected to the hemi- 
sphere of radius Sz/2 from the point (0,0, —Sz/2)^ and from 
there projected to the plane kx — —Sz/2 from the point 
(s./2,0,0) 



giving the action 



27r A, 



' Pz{z) dz 
{s, - qz/2) z^ ■ 



(41) 



The endpoints of the integral are two of the positive roots 
of the cubic Psiz). We discuss the q > and ^ < cases 
separately. 



1. q<0 

For q < the level sets of the action are shown in 
For < 0, we see that Is is not smooth along 
0. To understand the origin of this behavior, let 
us consider the roots of Ps{z) as Sz 0. Two roots 
are positive, and one negative (as should be clear from 
Fig. [4|. Further, two of the roots are 0(5^), being the 
roots of the quadratic 



Fig. [10 

s. 



2hz^ -2sz{2^q)z^4sl, 

while the remaining root is Ah/szq^ and diverges to ±oo 
depending on the sign of h. For h > the two positive 
roots are those that are vanishing. Since these are the 
limits of the integral, we can see immediately that Is is 
0{s^). Further, the first correction is at order s^, so Is is 
smooth about Sz = 0. For < 0, one of the positive roots 
is diverging, corresponding to a trajectory that encircles 



FIG. 10. Level sets of the actions /3/iV (red) and Is = I3/N + 
Iszl/^ (blue) for q = —0.3. The first is cusped below /i = 0, 
the second above. 
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FIG. 11. Branch cuts (red segments) of the integrand and 
integration contour (dashed) in Eq. (41 ) for ^ < 0, /i < 0, and 
Sz ^ 0. z< and the unmarked root are both 0{sz). 



the tip of the cone. It is this divergence of the upper 
limit of the integral in Eq. ( [41] ) that is responsible for 
the cusp. To see this, note that the integrand has four 
square root singularities to be joined up by two branch 
cuts. For ^ < 0, there is one branch cut on either half 
of the real axis. The integral for the action is half the 
integral circulating the branch cut on the positive side 
(see Fig. 11). We can deform this contour so that it 
circulates the branch cut on the negative side and the 
pole at z = - both giving a contribution of 0{s^) - and 
the pole at infinity. The contribution of this latter pole 
gives the cusp 



FIG. 12. 
q = 0.3. 

h: 



\si + q. 



y 



Level sets of h on the half-plane for Sz — 0.45, 
The dark line is the separatrix corresponding to 



2. q>0 



For q > we look for interesting behavior associated 
with the separatrix h = \sl ^ q. In the half plane 
model the phase space corresponding to the truncated 
hyperboloid (/cq < 1/2) is a disc bounded by + (z — 
2/5^)^ = 4(s~^ — 1) (circles project to circles), and the 
separatrix consists of this circle plus the chord z = 2sz/q 
(see Fig. 12), corresponding to the square root divergence 
in the action integrand Eq. (41). 

The roots of Psiz) are all positive (in Fig. [2] the 
'third' solution lies on the other branch of the hyper- 
bolas of constant energy, outside of the physical phase 
space ko < 1/2). The endpoints 2;<^> are the greatest 
two roots at energies below the separatrix and the small- 
est two above it. These considerations show that the 
action is discontinuous at the separatrix. Repeating the 
analysis of the ^ < case shows that below the separatrix 
the action again develops a cusp 



h{H,S,)^h{H,0) 



2 



e 



■ q — h 



0{Sl). 
(43) 



2 



Q{-H) + 0{Sl), (42) 



where we now restore the modulus sign. Since dl^/dH < 
(see Fig. [4|, this is consistent with the contours in 
Fig. [To] The implications of the cusp will be discussed 
in Section IIIIDI Note that if we define a new action 
= Is + l'S'21/2, this has a cusp for > instead. 



D. Rotation angle and monodromy 

More significant than the action is its derivatives, 
which give us the period lattice vectors from Eq. (29). 
Specializing to the actions of Eq. ( 33 ) we find 



27r 



1 

1 

dis dl3 dis 

, ON dS^ dH , 



(44) 
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471- 

2ti- 





T 2Tt 

FIG. 13. (Left) As we circle the origin in Sz, H space for 
^ < the period lattice is deformed continuously, returning 
to its original form, but after shifting the lattice vector cor- 
responding to h by 27V in the direction. (Right) Schematic 
illustration of the rotation angle. While executing a single 
period T of motion on the reduced phase space the system 
rotates by an angle $ 



The vector 63 = 27r (^f^ |^ f^^ tells us how to exe- 
cute a closed orbit around the third circle of the three- 
torus: we evolve for a time 27r^ (this is then the period 
of the motion on the reduced phase space), change the 
overall phase of the spinor by 27r^, and rotate about 
the z-axis by 2ti^^. The rotation angle 



^{S,,H) = -2tt 



dh 



(45) 



is therefore the rotation about the z-axis associated with 
one period of the reduced motion (see Fig. 13). Com- 
paring with Eq. (42) we arrive at the surprising con- 



clusion that for g < 0, ^{Sz-,H) is not a single-valued 
function, but rather increases by 27r upon encircling the 
origin H = Sz =0. By contrast the period 



27r 



dh 
OH 



(46) 



is single- valued (though logarithmically diverging as we 
pass through the origin). T may be expressed as an el- 
liptic integral [18]. Note that for ^ > the separatrix 
divides the phase space into two disjoint regions (see 
Fig. 12). In each of these regions action- angle coordi- 



nates can be introduced without difficulty. 

The non-trivial mapping of the period lattice into it- 
self upon encircling the origin in (Sz^H) space is the 
characteristic signature of monodromy (see Fig. 13), and 



by Eq. (24) corresponds to angle variables that are not 



single- valued. The mapping of the period lattice vectors 
is written as 



^ e = Mi^ 



(47) 



where M is a integer- valued matrix of unit determinant 
(an element of the group SL{3, Z)) called the monodromy 
matrix. In our case 



M 




(48) 



Note that we can focus on the 2 — 3 subspace of the pe- 
riod lattice. Nothing interesting happens in the direction 
corresponding to /i = A/", reflecting the fact that N was 
scaled out of the problem. 

What is special about the point = = for ^ < 0? 
Recall from the discussion of Section [lII A| that for Sz = 
the reduced phase space is a cone and for < the tra- 
jectory encircles the tip, while for > it does not. 
From Eq. (31) the tip of the cone corresponds to the 
state (ai, ao, a_i) = (0, e~*^, 0), which is invariant under 
rotations about the z-axis. Thus the torus is pinched at 
this point: the circle corresponding to rotations through 
has contracted to nothing. Without such a singularity 
there would be no distinction between paths that circuit 
the origin and those that do not, and hence no possibility 
of non-trivial monodromy. Further, the structure of the 
singularity - known as a focus-focus singularity in the 
mathematical literature - actually fixes the monodromy 
without the need for explicit calculation of the actions 
[23H25j . To see this, let us consider the quadratic Hamil- 
tonian in the vicinity of the tip of the cone. After fixing 
ao = 1, Eq. ([5| gives 



^quad 



(49) 



In the range —2<q<0 this corresponds to an 'inverted' 
complex oscillator, as may be seen by defining 



1 f2_ 



" a/2 V 191 

i ( \q\ 
\/2 ^2-|q| 



1/4 



(ai +a*_^] 



1/4 



(50) 



with {2,tx7*}pg = 1. In terms of these variables 

/iquad = ^ (1^1' - \^?) (51) 



where O = \/|<7|(2 — \q\). The unstable and stable modes 
are then 



1 / A 



(52) 



satisfying {a+,al}pg = 1. In terms of these modes 



I PB 

Sz = o*_aj^ — a^a+. 



(53) 



The linearized equations of motion are a± = ±r^a±, 
showing that as a+ grows exponentially, a_ decays so 
as to conserve a^a_. We now recapitulate an argument 
from Ref. [26^ that shows how these simple considerations 
fix the monodromy. 

From Eq. (53), varying the overall phase of a^a_ 
amounts to circling the origin = = 0. When ei- 
ther a+ or a_ vanishes (and the other is small), we are 
on the pinched torus, the two components corresponding 
to the stable and unstable branches respectively. These 
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t = 



t = T 



( re 




FIG. 14. After one period of motion on the reduced phase 
space the system returns close to the pinch, but rotates by 
the rotation angle The angle between a+ and a- remains 
constant. We can circle the origin /i = = by winding 
a-(0) as shown, keeping a+(0) fixed. By the argument in the 
text, a+(T) must wind in the opposite direction, showing that 
the rotation angle changes by 27r. 



two branches are of course connected away from the lin- 



ear regime. If we start from |a-(0)| <C |tt+(0)| = r] 
at time t = for some <C 1 the system is close to 
the unstable branch of the pinched torus and will evolve 
< T to be close to the stable branch 
\a-{t^)\ = T]. The key observation is 
that the limit a_(0) ^ is well behaved: we approach 
the pinched torus but the evolution ^ t' excludes the 
'pinch'. Thus the phase of a-{t') will not wind with the 
phase of a_(0). However, since a^a_ is conserved, a+(t') 
must wind oppositely to the phase of a_(0). After evolv- 
ing for an additional time Vt~^ In | ^^Iq^ | to give a total 
time of one period T of the reduced motion, |a+| grows 
back to its initial value 77, and |a_| decays to its initial 
value, as the system passes close to the pinch. The phases 
of a± will have changed, however, and this change is just 
the rotation angle (as should be clear from the definitions 
in Eq. (50)). The winding of a+(T) in the opposite sense 



in finite time t' 
i.e. |a+(t')| <C 



to the winding of a_(0) corresponds to a 2tt change of 
the rotation angle, see Fig. [l4j 

We close this section by giving a more explicit illustra- 
tion of the rotation angle. If we fix the ao component of 
the spinor to be real, so that ao = ^/Y^^\arJ^^^^\aZl^ ^ 
the Gross-Pitaevskii equations become 



\ai\ 



[aia_ 



+ aia-i + |aip + |a-ip)a±i] + qa±i. 



(54) 



(we measure time in units of 62^). 
verse magnetization 5+ 

and large (approaching tt) rotation angles is shown in 



The trans- 
isy = V2{al + 
The evolution of 5+ for small 



Fig. 15 



E. Semiclassical quantization 

With this extensive groundwork laid we can finally dis- 
cuss the quantization of the problem. The semiclassical 
prescription of Einstein, Brillouin and Keller (EBK) is to 
quantize the actions according to 



li = (^i + Mi) ^7 Hi e Z 



(55) 



where the jUi are known as Maslov indices. In the present 
case there are no non- trivial Maslov indices, and we have 
the quantization rule (recalling that we set h = 1 in writ- 
ing down the actions in Section IIIB) 



/3 = — — /cosh6>#GZ^ 

27T 2 J 



(56) 



(the other two rules simply quantize N and Sz in the 
familiar way). The integral is related to the hyper- 
bolic area enclosed as explained earlier. The integer 
contours of the action are shown in Fig. [16) compared 
with the result of numerically diagonalizing the Hamilto- 
nian Eq. (|5|. A more accurate semiclassical quantization 



could be achieved using the results of Ref. [27 , but we 
do not require it here. 

Above the origin the eye naturally picks out negatively 
curved rows in the spectrum, and these correspond to the 
contours of the action 1^. Beneath the origin these rows 
are not so evident, reflecting the fact that the Is is not 
smooth here. Instead the eye picks out positively curved 
rows corresponding to the action = Is + \Sz\/2. 

An elegant description of the relation between the 
quantum spectrum in the semiclassical limit and the 
monodromy of the classical system has been given by 
Zhilinskii [28 . The EBK quantization rules tells us that 
locally the eigenstates form a lattice in the space of con- 
served quantities Fj j = 1^...N with lattice vectors 
given by changing each of the li by h 



dFj 



(57) 



Ik^Ii fixed 



(c.f. Eq. (28)). Thus the lattice of quantum states is just 
the reciprocal lattice of the period lattice. Accompany- 
ing the mapping of the period lattice vectors on circling 
the origin in Eq. (47) is the corresponding map on the 
'quantum' lattice 



(58) 
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FIG. 15. Evolution of transverse magnetization for q = —0.3. (Left) Initial conditions ai — 0.4, a_i = 0.3 (Right) Initial 
conditions ai = 0.4, a-i = —0.3, rotation angle approaching +7r 




FIG. 16. Spectrum of Eq. ^ with iV = 20 particles and q = 
—0.3. The dark lines are the integral contours of the action 
I3. Transporting an elementary cell of the lattice around the 
origin leads to non-trivial monodromy. 



with 



IV. SOLUTION OF THE QUANTUM 
HAMILTONIAN 

A. The Bethe ansatz equations 

In Ref. [12 , Bogoliubov used the algebraic Bethe 
ansatz (ABA) to solve the Hamiltonian Eq. ([5|, based 
on the mapping to hyperbolic spins [29|. An eigenstate 
of the Hamiltonian is written as 



Nr 



i*)=n 



K+ 



B+ 



\j \j + Nq/2 



0,iyK)K®\0,VB)B (60) 



where the Xj j = 1, . . . Nr satisfy the equations 



Nr 



A, + Nq/2 ^ A, - A, 



(61) 



Recall that = ^ (I'S'^I + 1) and vb = I, f according 
to the parity of N — Sz (see Eq. ([l2|). Since each factor 
in Eq. (60) creates a pair of particles the total number 



N is related to the number of roots Nr by 



N = 2Nr^\Sz\^2ub 



(62) 



'10 
{M-'Y = ( 1 -1 
,0 1 



The resulting 'defect' is illustrated in Fig. 16 



(59) 



These semiclassical considerations will be put on firmer 
footing in the next section, where we discuss the exact 
solution of the quantum problem. 



The energy of the state Eq. ( 60 ) is 



E{{\,}) = ^N{N-l)+q\S,\ 



4e 



Nr 



N 

i=i 



(63) 



(The more complicated form in Ref. [T2| is expressed in 
terms of Nr) 

In the next section we will give a derivation of these 
results that does not rely on the full machinery of the 
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ABA and is suited to semiclassical approximations. For 
the remainder of this section we discuss the solutions of 



Eq. (61) in quahtative terms. 



The simplest way to understand the character of the 
solutions is to interpret Eq. (61) as the extremal condi- 
tion of the 'potential' 



E 



InIA,; 



■UK\n\Xi+Nq/2\-\i]. 

(64) 

This corresponds to a set of Nr unit positive charges 
located at positions {A^} interacting among themselves 
a 2D Coulomb potential, with a pair of positive charges 
fixed at and —Nq/2 with strengths ub and uk respec- 
tively. Additionally there is a constant electric field push- 
ing the charges in a negative direction. Note that the po- 
sitions of the charges could be complex numbers, as in 
the related systems of equations considered in Ref.W, for 
instance. Since all charges in the present case are of the 
same sign, however, it is not hard to see that equilibrium 
configurations can only occur for all A^ real. 

Let's first consider the special case ^ = 0, which corre- 
sponds to zero magnetic field. In this case the spectrum is 
simple, consisting of SU (2) multiplets. This can be read- 
ily understood in terms of the above equations, where we 
now have a single fixed charge of strength i^b-^-^k at the 
origin. Let's start from a solution with given Sz- We can 
generate another with — 2, by adding one root (this 
keeps the number of particles fixed, see Eq. (62)) at the 
origin. The charge of the fixed charge decreases by one 
unit, but this is compensated by the new root, so the 
same configuration of the other roots is still a solution. 
Because the new root is at the origin, the total energy is 
unchanged (see Eq. (|63|). Proceeding in this way, and 
using both values of ub^ we can generate the whole SU (2) 
multiplet. 

It is straightforward to check that the correct ener- 
gies are reproduced in this case. Multiplying the Bethe 
ansatz equations Eq. (61 ) by Aj and summing over j gives 



(assuming that there are no charges at the origin) 



i=i 

giving the energy 



-N{N + 2vB + \S, 



(65) 



E 



62 {S,iS, + l)-2N, 



if N — Sz even 



2N \{S^ + l){S^+2)-2N, ifiV-5^odd 



(66) 



Which are the eigenvalues of ^ 



^ (S • S - 27V ) for total spin S = d^nd S = ^ 1 
I 



2N 



respectively. Note that Bose statistics limits S to even 
values for N even and odd values for N odd. 

At finite q the two fixed charges separate. When 
adding a root, reducing Sz by 2, one can choose to 
place it either between the fixed charges or in the re- 
gion A > max(— A^(7/2, 0) (some rearrangement of the 
other charges occurs). These two moves may be used to 
build up the spectrum and correspond to increments in 
two different choices for the action in semiclassical quan- 
tization (see Fig. 16). Moving between consecutive blue 
contours while staying on the same red contour corre- 
sponds to adding roots between the charges, and these 
go over to the SU{2) multiplets as ^ 0. Conversely 
moving between consecutive red contours while staying 
on the same blue contour corresponds to adding roots 
in the region A > max(— A/'(7/2, 0). In Section IV C 



we 



will verify that the cusps in these actions are reproduced 
correctly. 



B. Derivation of the Bethe ansatz equations 



Let us derive Eq. (61) without employing the full ma- 



chinery of the ABA. We begin by finding the discrete 
Schrodinger equation for the operator 

n = NqK^ + 2B^K^ - B^K' - B' . 

Writing an eigenstate as 



(67) 



n=0 



we obtain the following equation for the coefficients Cn 

- (TVi? - n + 1 + vb){Nr - n + VB)cn-i 

- (n + l)(n + 2vK)cn+i = (68) 

for a state with eigenvalue ^. Now we going to recast the 
problem as a differential equation for the polynomial 

^(A) = ^ Cn(-l)"A^^-^+"^ (A + 7V(7/2)^+"^ (69a) 
oc (A + Nq/2Y^ n " • (^^^) 

n 

With some lengthy algebra one can show that the dis- 
crete Eq. (68) is equivalent to the following differential 



equation for Eq. (69a) 



vk{vk-1) , vb{vb-1) 2{NR + VK + VB){\ + Nq/2) 



A2 



{X + Nq/2Y 



A(A + Nq/2) 



A(A + Nq/2) 



(70) 
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Now using Eq. (69b) we can check that if Eqs. (61) are satisfied we have a solution to Eq. (70) with eigenvalue 



(71) 



The relation between second order linear differential equations and equations of Bethe ansatz type is known as the 
Heine-Stieltjes problem. Finally we can use the Bethe equations again to show 



2NqiyK 



1 



A„ + Nq/2 



Nr{Nr - 1) 



Nr{vk + vb)- 



(72) 



After restoring all factors, plus the c-number pieces that we have dropped in passing to hyperbolic spins, we obtain 
the eigenenergy Eq. (63). 

For the purposes of semiclassical analysis, it is convenient to recast Eq. (70) as a conventional Schrodinger equation 
using the transformation x — so that x" ~ X — ~ 2^') e~^^ giving 



with 



x"+t(A)x = 0, 



t{X) = 1 + 



VBi'^B - 1) , '^KivK - 1) 



A2 



1 



1 



1 



(X + Nq/2)^ \X X + Nq/2 
1 



eigenvalue of B"+K° 

, ^ 

{Nr + i^B + Vk) 



Nq\X X + Nq/2 
P4X) 



{^-Nq{NR + VB + VK)l2). 



A2(A + A^g/2)2' 



(73) 



(74) 



where P4(A) is a fourth order polynomial. Notice that 
t{\) is just the eigenvalue of the transfer operator in- 
troduced in Ref. [12 . The resulting potential describes 
motion with combined 'centrifugal' and Coulomb poten- 
tials, the latter being of variable sign. 

The derivation of a Schrodinger equation for a related 
system was given in Ref. [31 using the method of sep- 
aration of variables, but we shall not elaborate on this 
connection here. 



C. WKB analysis of the Schrodinger equation 

An understanding of the spectrum for N large requires 
a semiclassical analysis of Eq. ( 73 ) by the WKB method, 
according to which the actions 



1 

2^ 



d\ V-t(A) = n,,5 + ^, (75) 



with Ua^h ^ When vb,k are not large the hanger 

modification I'B.Kij^B.K ~ 1) ^ {^b,k ~ Eq. (74) 



is required, as the naive WKB wavefunction behaves as 



X(A) 



(A + 7V(7/2)H[^^(^^-i)] 



1/2 



when A ^ 
when A -Nq/2 
(76) 



whereas the correct exponents are ub and uk respectively 
(see Eq. (|69|). 

The two actions 7^,5 correspond to motion between the 
two endpoints (roots of P4(A)) contained in the regions 



Ua : min(0, -Nq/2) < A < max(0, -Nq/2) 
7^5 : max(0, -Nq/2) < A < oo. 



(77) 



^^/—t{X) is the root density of the Bethe roots in the 
^ oo limit. It is possible to obtain the root density 
directly from the Bethe ansatz equations without using 
the Schrodinger equation (see Appendix [A]) . 

The above quantization conditions amount to fixing 
integer numbers of roots in the two regions. From 
Eq. (62) for the total number of roots, we have = 
—la -^\Sz\/2 + const., so that if la is smooth as ^ 0, 
then Ii) is not, and vice versa. 

To understand how the action can have a cusp, con- 
sider the q < case, in which case —Nq/2 separates 
1Za and 7^5. As 5'^ ^ the centrifugal potential at 
—Nq/2 is vanishing, leaving the 'Coulomb' part. For 
^ > (^ < 0) this Coulomb potential is repulsive (at- 
tractive) for A > —Nq/2 and attractive (repulsive) for 
A < —Nq/2. It is then not hard to show that the action 
corresponding to the region with the attractive potential 
has a cusp, as the turning point of the action 'falls in' to 
A = —Nq/2 diS Sz ^ 0. In the vicinity of this point the 



15 



integral Eq. (75) looks like 



5? 



NqX + Nq/2 



4(A + Nq/2)^ 

= -/ 



du\ 1 



SI 

n,2 



(78) 



where u 



^^^^Jv^^^^^ (note that we are concerned with 
Sz of 0{N)). This is appropriate to the side on which the 
turning point is falling in, because then the contributions 
from the other parts of the potential can be ignored. To 
find the contribution to the integral consider the deriva- 
tive 



slow dynamics of the condensate spinor, but a detailed 
theoretical description is presently lacking. Moving to 
still smaller systems could eliminate this complication. 

Turning to the quantum dynamics of a single mode, 
we note that the level spacing is set by the spin interac- 
tion energy 62- Typical magnitudes are 4.3 Hz for ^''Rb 
[1 and 33 Hz for ^^Na [3]. A direct observation of the 
quantized spectrum on these tiny energy scales seems un- 
likely at present. An alternative strategy is to ask how 
the dynamics is affected by this discreteness, leading to 
deviations from the classical predictions. The question 
was addressed in several recent papers that treated re- 
lated models [26l l32l [33j . 

The author acknowledges the support of the NSF un- 
der grant DMR-0846788. 
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Thus the corresponding action behaves as —\Sz\/2 as 
6'^ ^ 0: la has a cusp ioi E > {) and is cusped for 
E < {). The discussion at the end of Section |IV A| can 
then be sharpened with the identification (up to con- 
stants) Is = la = —I3 = —^3 — I 'S'^ 1/2. In this way we 
recover the monodromy discussed in Section [TH D 



Appendix A: Solution of the continuum Bethe 
equations 



The continuum limit of the Bethe Eqs. ([61| gives the 
integral equation for the root density p{X) 



= 



A A + Nq/2 

AG{A:p(A)7^0} (Al) 



V. CONCLUSION 

In this paper we have given a detailed analysis of the 
semiclassical dynamics and spectrum of a spin 1 Bose 
microcondensate in the single mode approximation. This 
simple system proves to be rather rich, displaying Hamil- 
tonian monodromy in the g < case and a separatrix 
that divides the phase space for (J > 0. Both of these 
classical phenomena have distinctive quantum analogues. 

It remains to make a few comments about the relation 
to recent experiments in ultracold gases. As mentioned 
in the introduction, several experiments have observed 
(semi-) classical dynamics consistent to some degree with 
the single mode approximation [1- 3 . Ref. [3 is of partic- 
ular interest for its use of Faraday rotation spectroscopy 
in addition to the usual Stern-Gerlach separation of the 
different spin components. While the latter is only sen- 
sitive to the relative occupancies of the m = +1,0,-1 
states, the former is capable of measuring the transverse 
magnetization, which in principal allows the rotation an- 
gle to be extracted (c.f. Fig. 15). All of these experi- 
ments observe quite significant damping of single mode 
dynamics, indicating that other modes may be signifi- 
cant. That the system size be small compared to the 
spin healing length (equivalently, the level spacing in the 
trap is large compared to the spin interaction energy) 
is a necessary but probably insufficient criterion for the 
validity of the single mode approximation. The temper- 
ature is typically much larger than the level spacing, so 
that many modes are occupied. A reasonable expecta- 
tion is that these 'fast' modes adiabatically follow the 



Write the root density as the jump in an analytic function 
p(A) = ^[/(A + iO)-/(A-iO)], 



SO that the integral in Eq. (Al) can be viewed as circling 



the branch cuts of the /(A). If /(A) has the form 



/(A) 



A(A + Nq/2) ' 



(A2) 



in terms of a fourth order polynomial Pa{\) (required 
to get the two regions of nonzero root density) one can 
evaluate the integral in terms of the residues at 0, —Nq/2 
and 00, which solves the problem if 



^/pm = 

VP4{-Nq/2) = 



2 

vxNq 



(A3) 



V'P4(A) A^, as A ^ cx). 



These conditions, together with a specification of the to- 
tal number of roots and the energy 



Nr 



j dXp{X) 



(A4a) 



E=^N{N-l)^q\Sz\-^ JdXXpiX), (A4b) 



2N 



fix the form of P4{X) uniquely. By evaluating Eq. (A4) 



in the same way one can show that the polynomial 



16 



arising from Eq. (74) solves the problem (in the large tical to the WKB conditions Eq. (75). The strength of 



N limit The only thing that doesn^t come out is 
the ^Nji{Nji — 1) in Eq. ([72|)). A more useful approach 
to generating the spectrum is to separately quantize the 
number of roots in each of the two regions where the root 
density is nonvanishing, and then evaluating the resulting 
energy. The two quantization conditions are then iden- 



the approach based on the Schrodinger equation is that 
the polynomial is explicitly given by Eq. ( 74 ) without the 
need to solve Eq. (A4). 



The observant reader may note that the form Eq. ( 74 ) 



is only consistent with the first two of the conditions 
Eq. (A3) in the limit of iyB,K ^ 1- The remedy is the 
Langer modification mentioned after Eq. ( 75 ) . 
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